Wobbling excitations in ^^^Dy and ^^^Yb 
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We study in the cranked Nilsson plus random phase approximation low-lying quadrupole excita- 
tions of positive parity and negative signature in ^^'^Dy and ^^^Yb at high spins. Special attention is 
paid to a consistent description of wobbling excitations and their identification among excited states. 
A good agreement between available experimental data and the results of calculations is obtained. 
We found that in ^^^Dy the lowest odd spin gamma- vibrational states transform to the wobbling 
excitations after the backbending, associated with the transition from axially-symmetric to nonaxial 
shapes. Similar results are predicted for ^^^Yb. The analysis of electromagnetic transitions, related 
to the wobbling excitations, determines uniquely the sign of the 7-deformation in ^^''Dy and ^^^Yb 
after the transition point. 

PACS numbers: 21.10.Re,21.60. Jz,27.70.-f q 



I. INTRODUCTION 



Deformation is an important ingredient of nuclear dy- 
namics at low energies [ll, Regular rotational bands, 
identified in spectroscopic data, are most evident and 
prominent manifestations of an anisotropy of a spatial 
nuclear density distribution. While an axial deformation 
of a nuclear potential is well established, there is a long 
standing debate on existence of the triaxial deformation. 
A full understanding of this degree of freedom in nuclei 
may give impact for other mesoscopic systems as well. 
In particular, the importance of nonaxiality is discussed 
recently for metallic clusters Q and atomic condensates 
In nuclei, the nonaxial deformation involves the ori- 
entation degree of freedom that could be exhibited as 
wobbling excitations P, 0, 0] and chiral rotation \§\ . 

Various models predict that nuclear ground states in 
the mass region 150 < A < 180 are characterised by 
axially deformed shapes. With increase of an angular 
momentum, numerous calculations demonstrate that de- 
formed nuclei undergo a shape transition from axially 
deformed shape to nonaxial one (cf Ref . Such a tran- 
sition may be manifested as a backbending, known since 
early seventies [3]- There are different mechanisms re- 
sponsible for the backbending, in particular, in the tran- 
sitional nuclei with N ~ 90 (cf RefsES, lll|, IjJ). For 
example, we found (see RefjlJ) that in ^^^Dy such a 
transition can be explained as a result of vanishing of 
the 7-vibrational excitations in the rotating frame, while 
in -'^^^Yb it is due to a rotational alignment of a neu- 
tron two-quasiparticlc configuration. In the both cases, 
neither mechanisms do not provide, however, a non am- 
biguous proof of the onset of triaxial shapes. 

Nowadays, according to a general notion, the analy- 
sis of specific low-lying excited states near the yrast line 
could shed light on existence of the nonaxiality. For non- 
axial shapes one expects the appearance of low-lying vi- 
brational states, that may be associated with a classical 



wobbling motion. Such excitations (called wobbling ex- 
citations) were suggested first by Bohr and Mottelson in 
rotating even-even nuclei [l| , and analysed soon within a 
simplified microscopic approach in Refsfl^. [isl (see also 
Ref, 16 and references therein). According to the micro- 
scopic approach [sl, Q , the wobbling excitations are vi- 
brational states of the negative signature built on the 
positive signature yrast (vacuum) state. Their character- 
istic feature is collective E2 transitions with AI — ztlTi 
between these and yrast states. First experimental evi- 
dence of such states in odd Lu nuclei was reported only 
recently Q- 

Properties of the wobbling excitations at different an- 
gular momenta may be studied within the asymmetric ro- 
tor model (ARM) IJ (see also Appendix) . The extension 
of this model for odd nuclei has been used recently for 
the analysis [l3] of experimental properties of the second 
triaxial superdeformed band in -"^^^Lu, that carries a few 
features associated with the wobbling excitations. The 
classical dispersion equation for the wobbling mode [l|, 
with irrotational moments of inertia, was used to describe 
a spectrum. The moments of inertia were fitted in order 
to reproduce the experimental data. As a result, a con- 
troversial 7-reversed moment of inertia was introduced 
in Ref[l3 to resolve difficulties with the data interpre- 
tation. To overcome this problem the alignment of odd 
quasiparticle has been taken into account later in Ref .lSl. 
It appears that the approach suggested in RefsflTl. Il8l 
may explain some tendencies, being, however, a crude 
approximation to a full physical picture of the observed 
phenomenon. 

To explain the same data in ^^^Lu, a non self-consistent 
microscopic analysis, based on the cranked Nilsson po- 
tential, was performed in RefsHol. [20I As a basic tool, 
the dispersion equation for wobbling excitations, derived 
in the time-dependent Hartree-Bogoliubov approach in 
Refj^, has been used. Based on the solution of the mi- 
croscopic equation, in Refill it was concluded that the 
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pairing correlations do not affect the wobbling excita- 
tions, and this should be considered as a specific feature 
related to this mode. On the other hand, it was also 
found that the wobbling excitations are very sensitive to 
a single-particle alignment. As it is known, the align- 
ment decreases the pairing correlations. Therefore, the 
question arises about the validity of this conclusion. Fur- 
thermore, the authors admitted that the kinematic mo- 
ment of inertia was not described properly (see the 
discussion in Refl20h. In addition, the analysis [191 [20| 
is based on constant mean field parameters that are not 
related to the physical minima of the chosen potential 
at different rotational frequencies. In fact, the deforma- 
tion parameters have been taken semi-quantitatively to 
reproduce roughly experimental data. This may be not 
so crucial for the description of strongly deformed nu- 
clei, when the mean field parameters evolve slowly with 
a variation of the rotational frequency. However, this ap- 
proximation becomes questionable for the analysis of the 
rotational bands in the backbending region, which are 
the main subject of our paper. Evidently, mean field pa- 
rameters can change drastically at the backbending (see, 
for example, Ref[l3f). 

We recall that wobbling excitations depend on all three 
moments of inertia that characterize the nonaxial shape. 
Therefore, a self-consistent microscopic description of 
the moments of inertia is a necessary requirement to 
the theory that pretends to provide a solid analysis of 
experimental data. Recently, we developed a practical 
method based on the cranked Nilsson potential with sep- 
arable residual interactions for the analysis of the low- 
lying excitations near the yrast line. In contrast to pre- 
vious studies of low-lying excitations at high spins (cf 
Refs.flEI, [TgI ITqI . [20j'). we paid a special attention to the 
self-consistency between the mean field results and a de- 
scription of low- lying excitations in the random phase ap- 
proximation (RPA). Hereafter, we call our approach the 
CRPA. All the details of this approach are thorou ghly 
discussed in ReflTsl. The analysis of Ml-excitations [2ll |. 
shape-phase transitions and the behaviour of the positive 
signature excitations at the backbending [IS] confirmed 
the importance of the self-consistency for the description 
of moments of inertia. 

In the present paper we will study the low-lying neg- 
ative signature quadrupole vibrations and, in particu- 
lar, the wobbling excitations in ^^^Dy and ^^^Yb that 
are subjected to the backbending at high spins. Note 
that a microscopic formulation of the wobbling exci- 
tations in the form of the macroscopic rotor model 
was done by Marshalek in the time-dependent Hartree- 
Bogoliubov approach for the spherically symmetric, 
pairing-|-quadrupole-quadrupole Hamiltonian in the prin- 
cipal axes (PA) frame [6,]. Following the same approach 
in the PA frame, Shimizu and Matsuzaki [l^ analysed 
the electromagnetic transitions and presented a few nu- 
merical examples, using the cranked modified harmonic 
oscillator (standard Nilsson) potential. In the standard 
Nillson potential, however, the inertial properties are not 



described correctly (see discussion in Reflish . From point 
of view of the self-consistency, in order to study the vi- 
brational excitations with this potential, it is appropriate 
to use the double-stretched quadrupole interaction ^2^. 
We will provide a refined microscopic description of the 
wobbling excitations for doubly-stretched quadrupole in- 
teraction and their electromagnetic properties in a time- 
independent, uniformly rotating (UR) frame. In nu- 
merical analysis we will pay an attention to the self- 
consistency between the mean field, vibrational excita- 
tions and their electromagnetic properties. We recall that 
in the UR frame one considers all negative signature exci- 
tations, including the wobbling one. To identify the wob- 
bling excitations in experimental data, a few criteria are 
proposed in literature, such as: a large collectivity, zig- 
zag behavior of the B(E2) transition probability from a 
given band into the yrast one with AI = ±1, when one of 
the transitions is almost dominant (23| . We will present 
a microscopic procedure that provides a definite answer 
how to identify the wobbling excitations. This procedure 
includes the analysis of inertial properties, B(E2)- and 
B(Ml)-transitions probabilities. 

The paper is organized as follows: in Section II we 
briefiy review the main details of our approach that is 
thoroughly discussed in Reffisl. In Sec. Ill we study the 
lowest negative signature RPA excitations. The main fo- 
cus of this section is devoted to the definition of the spe- 
cific characteristics, associated with the wobbling mode. 
This task will be done on the basis of the comparison 
the CRPA approach and the ARM. The conclusions are 
finally drawn in Sec. IV. To complete the analysis, we 
review properties of the wobbling mode in the ARM in 
Appendix. 



II. THE MODEL 



A. Basic properties of the mean field 

Our description is based on the Hamiltonian defined in 
the UR frame 

Hn ^ H - hnX ^ Ho- KNr - TiVlJ^ + V (1) 

r—n,p 

The unperturbed Hamiltonian ~ '^ii^Nuii) + 

hadd{i)) consists of the Nilsson Hamiltonian 

T)2 1 

7 -T /22 22 22\ 

hNii = t: \- -m(U!.x-. + UJ2X2 + uj.,x^) 

2m 2 

-2K?iwool • s - K/x?iwoo(l^ - (l^)Ar), (2) 
and the additional correction term [241 
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The correction term restores the local Galilean invariance 
broken in the rotating coordinate system and improves 
the description of the inertial properties in the Nilsson 
model (see Reffl3l). The chemical potentials At (t =n or 
p) are determined so as to give correct average particle 
numbers {Nt)- Hereafter, (...) means the averaging over 
the mean field vacuum (yrast) state at a given rotational 
frequency fl. 

The two-body potential in Eci.((T]) includes the 
monopole pairing, doubly-stretched quadrupole- 
quadrupole and monopole-monopole interaction and 
spin-spin interaction. The Hamiltonian ^ possesses 
an inversion and signature symmetries. Using the 
generalized Bogoliubov transformation for quasiparticles 
and the variational principle (see details in Refl25l). 
we obtain the Hartree-Bogoliubov (HB) equations for 
the positive signature quasiparticle energies e,; (protons 
or neutrons). The positive signature (r — +1) state 
is defined according to the Bogoliubov transformation 



= T.A^k^ct +VuCk) and e-^"^a+e— ^'^ = -ia+ . 



Here, |fc) denotes a s.p. state of a Goodman spherical 
basis (see Refl26l). By diagonalizing the Hamiltonian 
at the rotational frequency fi, we obtain quasiparticle 
states with a good parity tt and signature r. It is enough 
to solve the HB equations for the positive signature, 
since the negative signature eigenvalues and eigenvectors 
are obtained from the positive ones 



denotes the negative signature (r 



(4) 



The index i 

state {e^'^^^ e~'^'^^'' = iaf). For a given value of the 
rotational frequency Q the quasiparticle (HB) vacuum 
state is defined as a^l) = aj|) = 0. 

We solved a system nonlinear HB equations for ^^^Dy 
and ^^^Yb on the mesh of deformation parameters /3 and 
7 defined by dint of the oscillator frequencies in Eq.([21) 



prescription [2 



A.(0) [l-kilh?] 
A.(0)i(%)2 



i}> fir 



(6) 



where flc is the critical rotational frequency of the first 
band crossing. 

It is well known that for a deformed harmonic oscilla- 
tor the quadrupolc fields in doubly-stretched coordinates 
fulfill the stability conditions 



/^ = 0,1,2 



(7) 



The tilde indicates that quadrupole fields are expressed 
in terms of doubly-stretched coordinates Xi — (uji/ujo) Xi 
and contain different combinations of the non-stretched 
quadrupole Qq oc {2z^ — x^ — y-^), Q2 oc — V3(x^ — y^) and 
monopole M oc operators quantized along the axis z 
(cf (2I1). The condition Eq.© holds, if the nuclear self- 
consistency 

u;!{xl) ^ u;l{xl) ^ u;i{xl) (8) 

is satisfied in addition to the volume conserving con- 
straint. In virtue of the condition ([8]) the doubly- 
stretched residual interaction does not contribute to the 
mean field results in the Hartree procedure. Enforcing 
the stability conditions Eq.© in the Hatree-Bogoliubov 
approximation, we search the HB minimum for the 
Hamiltonian (IT|) at a given rotational frequency. While 
the mean field values of the quadrupole operators Qq, Q2 
are nonzero, the doubly-stretched quadrupole moments 
(Qo) and (Q2) are zeros (see FigllJ for equilibrium defor- 
mations (see Fig [5]). The deviation from the equilibrium 
values of the deformation parameters /3 and 7 results into 
a higher HB energy, indeed. 
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The Nilsson-Strutinsky analysis of experimental data on 
high spins in ^^^Dy [231 indicates that the positive parity 
yrast sequence undergoes a transition from the prolate (a 
collective rotation around the z-axes) towards an oblate 
(a non-collective rotation around the x-axes). In order 
to compare our results with available experimental data 
on excited states (27i |. in comparison with our previous 
work [1^ , we extended the range of the values for 7 from 
7 = 60° (an oblate rotation around the y-axes) to 7 = 
—60° (an oblate rotation around the x-axes). At each 
rotational frequency, in each mesh point we calculate self- 
consistently the total mean field energy Ehb = (Hn). 
In the vicinity of the backbending, the solution becomes 
highly unstable. In order to avoid unwanted singularities 
for certain values of Q, we followed the phenomenological 
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FIG. 1: (Color online) The rotational behaviour of the cal- 
culated monopole and quadrupole moments. The "double- 
stretched" and standard values are connected by dashed and 
solid line, respectively. 

The results of our calculations conform to the results 
of the Nilsson-|-Strutinsky shell correction method (com- 
pare our Fig[2] with Fig. 3c in Reff27l). although we ob- 
tain slightly different values for the equilibrium deforma- 
tions. In the analysis of Ref 27 the pairing correlations 
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are omitted. In addition, in the Nilsson+Strutinsky shell 
correction method the rigid body moment of inertia sim- 
ulates the inertial nuclear properties, which is different 
from the microscopic one, even at high spin region (see 
below). Moreover, the use of the Nilsson+Strutinsky re- 
sults destroys the self-consistency between the mean-field 
calculations and the RPA analysis. Therefore, to keep a 
self-consistency between the mean field and the RPA as 
much as possible, we use the recipe described above. 




FIG. 2: Equilibrium deformations in /3-7 plane as a func- 
tion of the angular momentum I = (Jj,) — 1/2 (in units of 
h). The equilibrium deformations for ^^®Dy provide the lower 
mean field energies in the region — 7r/3 < 7 < (filled cir- 
cles) in comparison with those (open squares) obtained in 
Reflla. The both branches of the equilibrium deformations 
are obtained enforcing the condition Eq.((7]). The maximal 
difi'erence between the minimal HB energies at the positive 
and negative equilibrium 7-values does not exceed ~ IMeV 
for I'^SDy. 

The triaxiality of the mean field sets in at the crit- 
ical rotational frequency Jlc which triggers the back- 
bending in the considered nuclei due to different mech- 
anisms. As shown in Figl^l we obtain the critical ro- 
tational frequencies, at which the first band crossing 
occurs, h^lc ~ 0.25MeV (lO^i ^ 12K) for i^zyb and 
~ 0.3MeV {l^h l%n) for is^Dy. The contribution 
of the additional term, Eq.Q, was crucial to achieve a 
good correspondence between the calculated and exper- 
imental values of the crossing frequency in each nucleus 
(see discussion in Reflish. 

In ^^^Dy we obtain that the 7-vibrational excitations 
[K = 2) of the positive signature tends to zero in the 
rotating frame at the transition point, in close agree- 
ment with experimental data. At the transition point 
there are two indistinguishable HB minima with different 
shapes: axially symmetric and strongly nonaxial. It is in- 
teresting to note that this behaviour is symmetric with 
respect to the sign of the 7-deformation, although the 
difference between the HB energy minima for 7 — ±20*^ 
is about ^ O.SMeV. The increase of the rotational fre- 
quency changes the axial shape to the nonaxial one with 
a negative 7-deformation (7 ^ —20°). The transition 
possesses all features of the shape-phase transition of the 



first order. One may expect the appearance of the wob- 
bling excitations near the transition point, since after a 
shape transition there is a strong nonaxiality in ^^^Dy. In 
contrast with ^^^Dy, in ^^^Yb the axially symmetric con- 
figuration is replaced by the two-quasiparticle one with 
a small negative 7-deformation. There, the backbending 
takes place due to the rotational alignment of a neutron 
*i3/2 quasiparticle pair. The nonaxiality evolves quite 
smoothly, exhibiting the main features of the shape-phase 
transition of the second order. The question arises how is 
reliable our description or how self-consistently are done 
our mean field calculations ? 
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FIG. 3: (Color online) The rotational dependence of the 
dynamic S?^^ = —cPEhb/cI^^ = d{jx)/dQ. (a solid line) and 
the Thouless-Valatin Sjtv moments of inertia (a dash line). 
The experimental values SJ'-'^^ = 4/A_B^ are denoted by filled 
square connected by a thin line to guide eye. Here, l\E~, is 
the difi'erence between two consecutive 7-transitions, and E^, 
is the 7-transition energy between two neighboring states that 
differ on two units of the angular momentum 

One of the conclusive tests of the self-consistency mi- 
croscopic calculations is the comparison of the dynamic 

(2) 

moment of inertia S^^g , calculated in the mean field ap- 
proximation, and the Thouless-Valatin moment of inertia 
Stv , calculated in the RPA. They must coincide (see the 
results for exactly solvable model in Ref '29' ) , if one found 
a self-consistent mean field minimum and spurious solu- 
tions are separated from the physical ones. Our results 
(see Figl3]) demonstrate a good self-consistency between 
the mean field and the CRPA calculations, indeed (see 
also the discussion in Reffisl) . We emphasize that the in- 
clusion of the correction term, Eq.Q, is crucial to achieve 
a good description of the inertial nuclear properties. 



B. Negative signature excitations 



To describe quantum oscillations around mean field 
solutions the boson-like operators h'^^ = a^a^,b'^i = 
a^a^jbj^j = a^a^ are used. The first equality in- 
troduces the positive signature boson, while the other 
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two determine the negative signature ones. These two- 
quasiparticle operators are treated in the quasi-boson ap- 
proximation (QBA) as an elementary bosons, i.e., all 
commutators between them are approximated by their 
expectation values with the uncorrelatcd HB vacuum 
[30]. The corresponding commutation relations can be 
found in Refl25l. In this approximation the positive and 
negative signature boson spaces are not mixed, since the 
corresponding operators commute and — Iiu(r = 
+1) + li^(r = — 1). The analysis of the positive signa- 
ture term Ho^ir = +1) is done in Refill. 

In the UR frame, the negative signature RPA Hamil- 
tonian has the form 



A'3 = l,2 



(-)2 
Ms ■ 



(9) 



where = -|- {E-q = + e-^ are two-quasiparticle 
energies. Hereafter, we use the following definitions: the 
index \i runs over ij, ij and the index /is is a projection on 
the quantization axis z. The double stretched quadrupole 

operators q'\^ = CQi""* (C = w^^^/^o)' '§2"'' = ^'02~'' 
(?7 = uj^LOyjuj^) are defined by means of the quadrupole 

operators Q^m (ni=0,l,2) 



•2+rn + (r+3)/2 , 



^2(1 + 5^0) 



(10) 



where Q\m = f^Yxm- We recall that the residual doubly- 
stretched interaction does not distort the mean field de- 
formations found self-consistently for the Hamiltonian 
([1]). The violation of the self-consistency leads, however, 
to unreliable description of the vibrational states. 

The RPA Hamiltonian ([9]) contains only the isoscalar 
part of the quadrupole interaction, since we would 
like to establish a connection between the microscopic 
approach and the phenomenological ARM to see ev- 
idently similarities and differences. Moreover, in the 
considered nuclei the main contribution of the isovec- 
tor quadrupole-quadrupole interaction is located in en- 
ergy region around ~ 3MeV and is responsible for 
Ml excitations [2l[. One should keep in mind, how- 
ever, that the isovector part of the quadrupole inter- 
action may be important for the analysis of wobbling 
excitations in odd-odd nuclei, when there is a differ- 
ent orientation of the neutron and proton s.p. high-j 
orbitals and strong A/ 1-transitions are observed along 
yrast and/or yrare states. The spin-spin interaction, 
V„„ = -5ET=o,i«'-(^)Er=±Ep=o4(«M[r])^ (sec the 
discussion in RefsIT^. [21I ) is not essential for low-lying 
quadrupole vibrations and is omitted in our analysis. 
The pairing interaction does not contribute to the bo- 
son Hamiltonian, H^[r = — 1], since it is of the positive 
signature. Notice, however, that the matrix elements of 
the operators depend on the pairing interaction which 
affects the RPA solutions. 

The linear boson part of the double stretched operators 



has the following form 

Q['^ = -JE^(&^ + ^)' (11) 



Q2~^ = -iE/2M(&^-^)' /2M=W2M (12) 



Here, gi^, ^^e real matrix elements of the operators 

Qi~\ '52~^ respectively (the properties of matrix ele- 
ments of the operators involved in the Hamiltonian ^ 
can be found in Ref'2(?). 

We solve the RPA equations of motion 

[X,,P,,] = tS,,,, (13) 

where 

x, = ^x;(fe+ + M, F.-*^p;(6+-M (14) 

are, respectively, the collective coordinates and their 
conjugate momenta (hereafter, we use in all equations 
h ^ 1). The RPA eigenfunction 

1 



Ot\RPA) = ( A, - iP, ) \RPA) 



.,^h,)\RPA) 



(15) 



defines the amplitudes 1/'^'''' and ip^/i^^ by means of the 
generalized coordinate and momentum amplitudes. The 
ket vector | RPA) denotes the RPA vacuum (yrast state) 
at the rotational frequency Q. 

The solution of the above equations (|13p determines 
the generalized coordinate and momentum amplitudes 



with unknown coefficients 





f Xi?2 


Efif2ti 




El-0^ 




fXi?2 






El-0. 



(16) 



^1 =EAm^m =-71 ^.,Qi 
Ri - 2^ h^x^ = 



i(-) 



(17) 



and eigenvalues uj^. To find the eigenvalues 10^, one trans- 
forms the system Eg. lfTB)) to the form 



£'ii(w^) - - 
X 



R'(Di2{uJ,) + R2 



+ R^DuM - 



X 



= 



(18) 



The condition 

F(uj^) = det (D - -) = (19) 

X 

determines all negative signature RPA solutions. The 
matrix elements Dum{uJv) = fk4Lfni,p.C'^"'/{El-ujl) 
involve the coefficients C^™ — oj^ ior k ^ m and 
otherwise. Although the determinant has a dimension 
n = 2, one obtains a huge family of RPA solutions with 
different degree of the collectivity. Among collective so- 
lutions there are solutions that correspond to the shape 
fluctuations of the system. Notice that the direction of 
the angular momentum is fixed in the UR frame. 

III. THE WOBBLING MODE 

A. Janssen-Mikhailov equation 

We recall that the doubly-stretched residual interac- 
tion restores in the RPA the rotational symmetry bro- 
ken at the mean field approximation. Therefore, for the 
cranking Hamiltonian the conservation laws imply in the 
RPA order 

[Hn,JyTiJz] = ±^{Jy T iJz) (20) 

This condition is equivalent to the condition of the exis- 
tence of the negative signature solution ~ Vl created 
by the operator (see Reflsih 

Y^^-ll±£^^ f = (ft)t, [f,ft] = l (21) 
V2(Jx) 

This operator describes a collective rotational mode in 
the subspace of Hq{~) 

[ffo(-),ft]=f}ft, (22) 

arising from the symmetries broken by the external rota- 
tional field (the cranking term). 

This is true for a pure harmonic oscillator model only 
[2^. However, our Hamiltonian ([1]) contains the addi- 
tional term hadd, Eq.Q. Since the additional term hadd 
has a term proportional to Zj;-operator, the conservation 
laws ((20|) are broken for the Hamiltonian ([T]). Never- 
theless, these laws can be fulfilled in the RPA order by 
changing the strength constant in Eq. (fT9)) in order to ob- 
tain the solution = (IS] . To check this fact, we 
calculated the RPA secular equation (fT9|) for the mode 
uji, = 11, with and without the additional term hgdd, for 
the same strength constant (see also Fig. 10 in Ref ll3f) ) . 
The results evidently demonstrate that the conservation 
laws (PH)) . (1^ are fulfilled with a good accuracy (see 
Fig|4]). In fact, due to a smallness of the term propor- 
tional to the operator in the additional term, the vio- 
lation is almost negligible. 
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FIG. 4: (Color online) The evolution of the negative signa- 
ture RPA solution lj^ = 57 with (dashed line) and without 
(solid line) the additional term Eq.© as function of the rota- 
tional frequency, calculated at the equilibrium deformations 
(see Fig[2]). The straight dashed line, parallel to the h^l-axis, 
denotes the yrast line. 

Based on this fact, we exploit the conservation laws 
for the Hamiltonian that yield the following 
equations 

fij^ + E^jy = xe^/iM (23) 

njy + E^j; = xriBhf. (24) 

We will see below that these equations are fulfilled with 
a good numerical accuracy, since the redundant mode is 
separated from the other RPA solutions by our procedure 
for the strength constant. The parameters A and B 

- ( [g("' , iJy ]) = ^ ~h^.Jl = m2 + a/3Qo) (25) 
= { [Q(-' , J, ]> = ^ - 2,7(Q2> (26) 

are obtained with the aid of the commutator [J^ ± 
iJy,Qxm] = v/A(A -f 1) - m{m ± l)QAm±i (see Ref,3?). 
The above relations between the matrix elements, 
Eqs. (P5)) . ([Ml) . are a key point for the analysis of the wob- 
bling excitations at nonzero 7-deformation, i.e., {Q2) ^ 
0. Moreover, in virtue of the definitions for the phonon 
operator (Eqs.(Il5l),(I16l)) and the operator f (Eq. (|2Tl) ). 
employing Egs. f^ . iPH) . and Eg. ([55]) below, one can 
show that 

[f,a] = [f,6t]^0 (27) 

! — 1 

Following the procedure described in Ref 5, with the 
aid of Eqs. ((23|) . ((24| . one arrives to the equation 

{ojI - n^)A{uj,) = (28) 
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The determinant 
equations 



corresponds to the system of 



These relations can be associated with the system of 
equations 



for unknowns 



^^<>xz' 1 



R2 



that does not contain anymore the solution uj,y 
introduced above the following notations 



where 



TV /z 



(29) 



(30) 



VL. We 





— ^ — 


Oy 






^xz 


— ^ — 

— Ox 




ojls/n 


(31) 













(32) 



(jJ^Sztt — fl'^xyb ~ 

ft'^XzO' — OJySyh — 

for unknowns a ~ fiSy and b 



of definitions, Eqs.([3T 
following form 



(36) 

With the aid 
we rearrange this system to the 



^{"^x - 3. - ^uS^)a - uu^i'^y + uj,S^)b = (37) 



In virtue of the relations, Eqs. ([35|) . between the un- 
knowns a and 6, we define the effective moments of inertia 
in the UR frame 



'^2^ — ^y + — — 3y + ris" ^ 

bz 
a by 



(38) 



and Q^a: = (Jx)/^ is the definition of the kinematic mo- 
ment of inertia. From the system (|29p . one obtains the 
relations between unknowns and 



xy 



i^lfSy 



(33) 



which are helpful for our analysis below. 

The condition A^i,) — leads to the Janssen- 
Mikhailov equation 



A(w^) 



]s/n][^x-^z~Lols/n] 



[3y+J7S'][5^+OS'] 



(34) 

that determines all vibrational modes of the negative sig- 
nature excluding the solution uj^ — Vl. We stress that 
the solution of this equation alone is meaningless. While 
Ea. ([M|) does not depend on the strength constant, the 
violation of the conditions ((23| . ((24|) . by arbitrary varia- 
tion of the 7-deformation or pairing gap destroys the link 
between the systems, Eq. p^ and Eg. ipS)) . As a result, 
the redundant mode can not be removed from Eq. (jl9p 
and one is not able to obtain Eg. ((34)) . Providing the 
wobbling solution, this equation has, however, a different 
form than the Bohr-Mottelson classical equation. Below 
we present a simple derivation of the microscopic analog 
of the latter one. 



B. Marshalek moments of inertia 

The relations Eg. ([55)) are equivalent to the following 
ones 



'^iSy ^^xy (J-^I/Sy 



TnSz UJj^Sz ^'^XZ 



(35) 



that depend on the RPA frequency. The definitions 
are the same as those obtained by Marshalek in the time- 
dependent Hartree-Bogoliubov approach but in the PA 
frame Q. The determinant for nonzero solutions of the 
system Eq.([37]) yields the nonlinear equation similar to 
the classical expression for the wobbling mode pj] 



U^u—w — ^ 




"^3 J 



(39) 



with the microscopically defined moments of inertia, 
g Eq. (|55)) was obtained first by Marshalek in the PA frame 
from the equations for the amplitudes of the angular fre- 
quency time oscillations. In the UR (time-independent) 
frame we obtain this equation with much less efforts as 
well. 

It is evident that for the rotation around the axis x the 
wobbling excitations with different collectivity could be 
found from Eg. ()55)) . if the condition 



•Jx > ^2 1 ^3 I'''' 



^2 1^3 ) 



(40) 



is fulfilled. One may expect that for the RPA solutions, 
different from the wobbling mode, this condition should 
not hold. Notice that Eg. (fT9| contains the solutions to 
Eg. ()39p but not vise versa, since the above constraint, 
Eg. ([^D)) . valid for the later case, is not required for the 
former one. 

We obtain quite a remarkable correspondence between 
the experimental and calculated values for the kinematic 
moment of inertia for the both nuclei (see top panels in 
Figs l5l6p . We trace also the evolution of the irrotational 
fluid moment of inertia (see RefISC 



27r 



2tt ^ 



(41) 
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FIG. 5: (Color online) ^''^Yb. Top panel: the kinematic 
moment of inertia = {Jx)/^ (solid line), a rigid body 
(dashed line) and a hydrodinamical moments of inertia (dash- 
dotted line) are compared with the experimental values (filled 
squares). Experimental values ^x = I are connected by 
dashed line to guide eyes {TiQ, = E-y/2). Bottom panel displays 
the rotational dependence of the kinematic moment of inertia 
(solid line), Marshalek moments of inertia SJj^ (dashed line) 
and Ssg^ (dash-dotted line) for the first RPA solution v = 1 
obtained from Eq. (fT^ . 
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FIG. 6: (Color online) ^^''Dy. Similar to FigH 



and a rigid body moment of inertia (i — I, 2, 3) 



^ ■ 



= ^mAi?2 I 1 _ y' A/3 cos(7 - y 



(42) 



as a function of the equilibrium deformations (see Fig|2). 
The irrotational fluid moment of inertia 5^*'"'"^ does not 



reproduce neither the rotational dependence nor the ab- 
solute values of the experimental one as a function of 
the rotational frequency. The rigid body values provide 
the asymptotic limit of fast rotation without pairing. Ev- 
idently, the difference between the rigid body and the cal- 
culated kinematic moments of inertia in the both nuclei 
decreases with the increase of the rotational frequency, 
although it remains visible at high spins. At very fast 
rotation hVl > 0.45MeV the pairing correlations are re- 
duced due to multiple alignments, and, therefore, the 
difference is moderated. 
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FIG. 7: (Color online) The rotational dependence of the RPA 
solutions, obtained with the aid of Eq. (|19|l (solid line) and 
Ea. (l39|l (dashed line). The RPA solutions, obtained with the 
aid of Eq.^, fulfil the condition Eq. (|40)) . 



The calculated values of the Marshalek moments of in- 
ertia, Eq. (|55|) . signal (see the bottom panels, Figs l5l6p 
that one should expect the appearance of the wobbling 
mode after a shape-phase transition in ^^^Yb and ^^^Dy 
for the first RPA solution found with the aid Eg. (fT9|) . 
As was stressed above, the separation of the redundant 
mode is an essential point for the RPA wobbling the- 
ory, that secures a reliable analysis of the RPA modes. 
Additionally, to be assure in the self-consistency of our 
RPA calculations, we compare the solutions, that may 
be associated with wobbling excitations, from different 
RPA equations dUl) and ([39]). We recall that Eq.^ 
depends on the strength constant x ^-nd contains differ- 
ent RPA solutions including the redundant mode, while 
this dependence is removed from Ea.(|39p. Evidently, if 
the redundant mode would not be removed from Eq. p9|) 
by our choise of the strength constants, the conditions 
Eas. ([^5)l . ([M)) should be broken. As a result, the con- 
sistency between Eq. (fT9)) and Ea. ((39|) would be broken 
as well and these equations would provide different solu- 
tions. A nice agreement between the roots of Eq. (fT9| and 
Eq. ([551) (see FigO confirms the vitality and the validity 
of our approach. 

Does the constraint (f40|) guarantee the existence of a 
collective wobbling mode solution ? Below we formulate 
criteria how to identify the collective wobbling excita- 
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tions. 



This result yields the following expression for the sum, 
Eq.®, 



C. Criteria for wobbling excitations 



we can define the unknown 
" = b/Sz (which determine 



By means of Eqs. dST]) , ((38 
variables ri~ 



a/S, 



the wobbling mode) such that 



Sy b 

Sz 



3f ) 



IW2 Sz 



Sy 




(43) 



Here, we used the fact that the dispersion equation for 
the wobbling mode, Eq. (|39p . can be expressed as = 
{X)VW^, where 

W2 = (l/5f - 1/5.), 



(l/5f 



l/3x) (44) 



With the aid of Eqs.(138l) and the definition, Eqs.(13Tj), it is 
easy to show that the ratio Sz^2^ / Sy^^^ = 1. It results 
in the following relation 



(45) 



Our next task is to find a formal expressions for the un- 
knowns in terms of the functions VF2,3 that depend 
on the microscopic moments of inertia. 
By means of the inverse transformation 




6t = V2Y,x;{dl - 6.) + p;(6t + d,) 



(46) 



we can express the operators Q^-] (^^^ Eqs. ([TT|) . (|12p . 
(|17p ) via phonon operators 



(-) 



(47) 
(48) 



Let us use the fact that the components of the 
quadrupole tensor commute, i.e., the condition 







(49) 



Here, we use the notation D^^q = F and D^^q = Oi, 
for other vibrational modes. Taking into account the 
definitions Eqs.HID,® (see also Eqs.(IlTD,(12n]),(IlSl)), we 
obtain the exact definitions for the unknowns rp2 asso- 
ciated with the redundant mode: 



1 



H2 



^AV2 



1 



ilB r]BV2 



{X 



(50) 



2\/(J.) 



1'2 



''1 ^2 



1 



(51) 



Suppose, that in the l.h.s of Eq. ([5T|) the sum, defined by 
all physical solutions excluding the wobbling one, is zero 
due to a mutual cancellation of different terms. Below, 
we will see that it is, indeed, the case. As a result, we 
obtain the equation for the unknowns ■ Resolving this 
equation by dint of Eq. (|45|) , we obtain 



1 



2\/(J.) 



W2 
W3 



1/4 



1 



2 J (J.) 



W2 



1/4 



(52) 
wob- 



These expressions are similar to those of the 
bling mode in the Bohr-Mottelson model (see Appendix, 
Eqs. (|A.14p ). although the quantities ^^2,3 are determined 
by the Marshalek moments of inertia and (Jx) ~ I +1/2 
(the factor 1/2 is arising due to the RPA contribution 
of the redundant mode, see discussion in Ref lSSi) . Note 
that similar expressions were obtained in Refllp in the 
PA frame (with the quantization condition (J^) = I and 
some additional phases). 

To identify the wobbling mode among the solutions of 
Eq. (fT9|l it is convenient to transform Eq. (j49|) to the form 



^ C^, = 0, C^= 4:{Jx] 



y—all 

From Eqs.jnni), dSS]) it follows that 

^U — Q. — 1; ^U — W = 1 



(53) 



(54) 



Thus, solving only the system of the RPA equations for 
the quadrupole operators, Eq. fTO)) . the condition Eg. fSl]) 
enables us to identify the redundant and the wobbling 
modes. 



D. Analysis of experimental data 

The experimental level sequences for all observed up- 
to-date rotational bands in ^^"^Yb and ^^^Dy are taken 
from ReflU, which is regularly updated. The rotational 
bands are numerated as rotational bands Bl, B2... in ac- 
cordance with [si^l . All rotational states are classified by 
the quantum number a which is equivalent to our signa- 
ture r. The positive signature states (r = +1) correspond 
to Q = 0, since the quantum number a leads to selec- 
tion rules for the total angular momentum I = a + 2n, 
n = 0, ±1, ±2 . . . (cf Ref|8|). In particular, in even-even 
nuclei the yrast band characterized by the positive sig- 
nature quantum number r ~ +1 {a = 0) consists of even 
spins only. The negative signature states (r = —1) corre- 
spond to a = 1 and are associated with odd spin states in 
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FIG. 8: (Color online) ^'^^Yb. Top left panel: the rotational 
dependence of the positive signature RPA solutions with even 
spins (tt = +,a — 0). Number in a circle denotes the RPA 
solution number : 1 is the first v — 1 RPA solution etc. 
Different symbols display the experimental data associated 
with B1,B2... bands (the band labels are taken in accordance 
with the definitions given in Ref lsih . Top right panel: the 
rotational dependence of the negative signature RPA solu- 
tions with odd spins {n — +,a = 1). The redundant mode 
u!,^ = Q is denoted as "0" and is displayed by the dotted line. 
Bottom panel: the rotational dependence of the coefficients 
^ R1R2 (see Ea. (|53p ) that are determined by the solutions 
of Eg.lfTOll. 



even-even nuclei. All considered bands are of the positive 
parity n — +. 

Before to make a comparison of our results with the 
data, let us analyse RPA solutions in order to identify 
a wobbling mode. Top right panels of FigsEl El display 
four lowest RPA solutions of Eq. (fTO)l and the redundant 
mode of Eq. p9)) as a function of the rotational frequency. 
We recall that these solutions are found at different equi- 
librium deformations (see Figl2]). Indeed, in the both 
nuclei the criteria Eq. (|54p uniquely determines the re- 
dundant and the wobbling modes. In FigslH El the re- 
dundant mode is manifested as a straight line (see the 
top right panels), while the corresponding coefficient cn 
is always -1 (the bottom panels). The rotational mode 
is separated clearly from the vibrational modes. Notice 
that the solutions, that are different from rotational and 
wobbling modes, contribute to the sum Eq. ([5T|l with a 
zero weight, as it was proposed above. 



FIG. 9: (Color online) ^^^^Dy. The same as for Fig|8] 



To compare our results with available experimental 
data on low lying excited states near the yrast line (s^l 
we construct the Routhian function for each rotational 
band v {v = yrast, [3, 7, ...) 



E,{I + 1)-E„{I-1) 
2 

(55) 



and define the experimental excitation energy in the ro- 
tating frame Tiw^^P = R^{VL) - Ryri^) [36]. In ^^^Yh 
it is known only one negative signature 7-vibrational 
state. The first RPA solution {1/ = 1) is a negative 
signature gamma-vibrational mode (with odd spins) till 
h^l w 0.28MeV. With the increase of the rotational 
frequency it is transformed to the wobbling mode at 
hfl « 0.32MeV (according to the criterion Eq.jSll)). The 
other solutions {ly — 2,3,4) contribute to the sum Eq. ([5T|) 
with a zero weight. Our results for v — 1 solution 
may be used as a guideline for possible experiments on 
identification of the wobbling excitations near the yrast 
line. Although the positive signature states have been 
discussed in Ref, 13, for completeness of our analysis we 
compare RPA results for the positive signature with up- 
dated database (34]. According to our analysis, the first 
RPA solution of the positive signature may be identified 
with /3-excitations at small rotation TiO < 0.2MeV. With 
the increase of the rotational frequency a strong mixing 
between 7- (the second RPA solution at low spins) and 
/3-excitations takes place. At hfl > 0.2 Mev the first 
RPA solution of the positive signature is determined by 
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a single two-quasiparticle neutron configuration (see dis- 
cussion and Table 1 in Rcf , 13). 

In ^^^Dy the first {1^ = 1) positive signature RPA so- 
lution carries a large portion of quasiparticle states with 
K = 2. Since the quantum number K is reliable at small 
angular momenta, we associate the first positive signa- 
ture RPA solution with 7-vibrational mode. Rotation 
leads to a strong mixing between the first and the sec- 
ond RPA solutions. With the increase of the rotational 
frequency the first solutions is separated, however, from 
the second one, while the latter one strongly interacts 
with the third RPA solution. At HQ « O.SMeV there 
is a crossing between Bl (yrast) band and B8 (excited) 
band which becomes the yrast after the transition point. 
The first positive signature RPA solution describes this 
transition with a good accuracy: the RPA mode goes to 
zero at Tiil w O.SMeV. We recall that according to our 
analysis [3| i namely this vibrational mode is responsible 
for the backbending phenomenon in this nuclei. 

The first negative signature RPA solution in ^^^Dy 
can be associated with the negative signature gamma- 
vibrational mode with odd spins. After the transition 
from the axial to nonaxial rotation, at h^l « 0.3MeV, 
according to the criteria Eq. (|54)) . the first negative sig- 
nature RPA solution describes the wobbling excitations. 
The mode holds own features with the increase of the ro- 
tational frequency up to fiil « 0.55MeV. There is a good 
agreement (see FiglS]) between Routhians for this band 
and for the experimental band BIO (or (-I-, l)i band ac- 
cording to Ref. UtI). On this basis we propose to con- 
sider the BIO band as the wobbling band in the range 
of values 0A5MeV < nn < 0.55 MeV {33n < I < 39n 
for this band). Note that the band BIO contains the 
states with 31 7i — 53 ?i. However, our conclusion is reli- 
able only for the states with / = 33?i — 39 ?i (or up to 
nn < 0.55 MeV). At hfl « 0.55 MeV a crossing of the 
negative parity and negative signature (positive simplex) 
B6 band with the yrast band B8 is observed. There- 
fore, for nn > 0.55 MeV (or for / > 39?i for the BIO 
band) one may expect an onset of octupole deformation 
in the yrast states. Since the octupole deformation is be- 
yond the scope of our model, based on the quadrupole 
deformed mean field, this feature will be discussed in 
forthcoming paper. 

The proposed criterion Eg. ([5^ is a necessary but not 
a sufficient condition to reach a conclusion that we have 
found a solution, related to the wobbling excitations. It 
is brought about by the formal equivalence between the 
classical Eq. (jA.3p and the microscopic Eq. ([M)) equations 
for the wobbling mode. Notice that our solution is de- 
termined in the UR frame, where the fluctuations of the 
angular momentum are absent (which arc responsible for 
the wobbling mode in the PA frame). To identify the 
wobbling mode we have to specify also the relation be- 
tween electromagnetic transitions in the Bohr-Mottelson 
model, defined in the PA systems, and our model, defined 
in the UR frame. 



E. Electromagnetic transitions 

Transition probabilities for XX transition llv 

> between two high-spin states are given by ex- 
pression 

B{XX-Jv ^ (56) 
{I IXlJ.i\I' <u'\M{XX■,^J.l -I)\u> |2 

At high spin limit (/ >> A, /' >> A), the transition 
from a one-phonon state into the yrast line state takes 
the form [37| 



B{XX;Iv^I'yr))^ 
\{RPA\ \m^^\XX-hi 



(57) 



\RPA)\ 



Here, M.^^'' {XX^ii) is the linear boson part of the corre- 
sponding transition operator of type X , multipolarity A 
and the projection /ii onto the rotation axis x in the UR 
frame. The commutator in (|57p can be easily expressed 

in terms of phonon amplitudes '4>']i'\ (/'IT'' (see Eg. (ITS]) ). 

With the aid of the transformation from x-to z-axis 
quantization (37| 

M{XXix,) = ^M3m(0' J'O) M{X\tX3), (58) 



and the definitions Eqs.(|T7 | . ([30|) . taking into account 

that {v\M2'^l=o.2\^) = (^^^^3=0,2) holds in the first RPA 
order, we obtain 



B{E2;Iv ^ I±lyr) - ( M^Ui^Oi ) 



V2 



V2 



(59) 



Here, we use M^^^ = {eZ/A)M. In virtue of Eq.dS!]), we 
arrive to the expression for the quadrupole transitions 
from the one-phonon wobbling state to the yrast states 



B{E2;Iw I±lyr) 



1 



4(Jx) 



W2 

1^3 



W2 



B 



(E) 



(60) 



which is similar to Eg. (jA.13p . Thus, we provide a com- 
plete microscopic definition of the wobbling excitations 
in accordance with the criteria suggested by Bohr and 
Mottelson for the rigid rotor [l|. Notice a clear differ- 
ence between the microscopic and the rigid rotor models: 
the microscopic moments of inertia, Egs. (j38p . should be 
calculated for the RPA solutions of Eg. (fTO|) that must 
obey to the condition Eq. 
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For the intraband transitions we have (see Eq.(43) in 
Refd) 



For the intraband transitions we obtain 



B{E2]Iv 
1 



I-2v) 



[v\M{E2-yi = 2)|i/> 



[E) 
2 



(61) 
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FIG. 10: (Color online) Reduced transition probabilities 
B{E2; I yr I — 2yr) along the yrast line. Experimental 
data (filled squares) are connected with a thin line to guide 
eyes. The results of calculations by means of Eq. (|66|) and of 
Ea. (|61|l are connected by dashed and solid lines, respectively. 

For illustrative purposes, to get a rough idea on the 
major trend of the quadrupole transitions, we employ 
the relations from the pairing-plus-quadrupole model (cf 
Refi) 



muiQp sin'-f 



xiQo) 
x{Q2) 



(62) 



between the deformation parameters /3 and 7 and two 
intrinsic quadrupole moments. By dint of Ea. (|62p and 
the definition of the quadrupole isoscalar strength x — 
47rma;g/5(r2) « iirmw^ / (SAR"^) {R « 1.2Aifm) one 
can express the £'2-transition probability via the defor- 
mation parameters /3 and 7 



I ±lyr) = e X 




(63) 



(64) 



where 8 = {9/lQTT^)e^Z^R'^. Depending on the sign of 
7-deformation (for W2.3 > 0) one obtains selection rules 
for the quadrupole transitions from the one-phonon wob- 
bling band to the yrast one 

a) i^or - 60° < 7 < : (65) 
B{E2;In^ I - lyr) > B{E2; I ^ I + lyr) 

b) For < 7 < 60° : 

B{E2;In^ ^ I + lyr) > B{E2; I ^ I ~ lyr) 



B{E2; In^ ^ / 



2n„)-ie/3W(J-7) (66) 



One observes from Eq. that for the transitions along 
the yrast line (n^, = 0) the onset of the positive (nega- 
tive) values of 7-deformation leads to the increase (de- 
crease) of the transition probability along the yrast line. 

Experimental values of B{E2, Iv I' yr) are de- 
duced from the half life of the yrast states f34| using the 
standard, long wave limit expressions B{E2,i ^ f) — 
Pii f) / (1.223 X 10^ E^) (e^ fm*) Here, the tran- 
sition energy E^ is in MeV and the absolute transition 
probability P{i f) = In2/T{i -> /) is the related to 
the half life r(i f) (in seconds). When comparing our 
results with experimental data, we take into account the 
Clebsh-Gordon coefficient (see Eq. ([55| ) up to / < 10?i. 
For / > 10?i the asymptotic value for the Clebsh-Gordon 
coefficient, which is 1, is used. Note that in the vicinity 
of the backbending the mean field description becomes 
less reliable (cf Reffssl). While the cranking approach 
should be complemented with a projection technique in 
the backbending region due large fluctuations of the an- 
gular momentum (cf Ref.'sO), its validity becomes much 
better at high spins. Evidently, the larger is the rota- 
tional frequency the stronger is a predictive power of the 
CRPA, since it is based on the cranking approach aimed 
for the high spin physics [Sj]. 

Experimental data for the quadrupole transitions along 
the yrast line are compared with the results of calcula- 
tions: (a) by means of Eq. ([CT|) and calculations (b) by 
means of Eq. (|M)l (see Fig ITU)) . In the calculations (a) we 
use the mean field values for the quadrupole operators. 
The calculations (a) evidently manifest the backbending 
effect obtained for the moments of inertia (see Figs l5l6p 
at hflc ~ 0.25, 0.3 MeV for ^^'^Yh and ^^^By, respec- 
tively. Thus, the use of the self-consistent expectation 
values {Qin ) is crucial to reproduce the experimental 
behaviour of the yrast band decay. The calculations (b) 
(Ea. (p5)l ) reproduce the experimental data with less ac- 
curacy, while providing the major trend of the transitions 
with the sign of 7-deformation. The agreement between 
calculated and experimental values of intraband B{E2) 
transitions along the yrast line is especially good after 
the transition point. 

At small rotational frequency, in the both nuclei, tran- 
sitions probabilities from the first positive and negative 
signature RPA solutions are much weaker compare with 
the quadrupole transitions along the yrast line (see FigfTUl 
and top panels of Figs lllll2p . At Ti^l - 0.05MeV the 
transition strengths for the first positive (r = -fl) and 
negative (r — —1) RPA solutions are : ~ 330e^/TO^ for 
i62Yb and - 500eVw* for ^^^T)y, with small differences 
between different transitions due to the Clebsh-Gordon 
coefficients. On the other hand, we obtain a good corre- 
spondence between the shape evolution and the selection 
rules for the both nuclei (see top panels of Figs fTlll2l 
and Figl2]). The transition probabilities, Eq. (l55)) . are cal- 
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FIG. 11: (Color online) ^^^yb. The electric B(E2)- (top) 
and the magnetic B(M1)- (bottom) reduced transition prob- 
abilities from the one-phonon bands to the yrast band. The 
positive (negative) signature phonon band is described by 
the first r = -|-l(r = — 1) RPA solution. One observes a 
strong dominance of the B(E2)- and B(Ml)-transitions from 
the wobbling states (r = — 1) with spin I to the yrast states 
with spin /' = / — 1, starting from the rotational frequency 
hi} > 0.28MeV. The transitions, calculated by means of the 
i/'p"^^' and ^Ir^^' phonon amplitudes, are connected by solid 
lines. On the right panels the results obtained by means of 
Eqs.dSOj, lis (with the aid of the variables 14^2,3, Ea. (|44)l ') 
are connected by thin lines, starting from the rotational fre- 
quency hfl ~ 0.3 MeV. This point is associated in our analysis 
with the appearance of wobbling excitations. 



culated by means of the ^p^i and (j))^' phonon amplitudes 
expressed via the coordinate and momentum amplitudes 
(see Eqs.(Il5l),([Tn]),(Il71),(ISni))- We compare these results 
for the first negative signature RPA solution (which is as- 
sociated with a wobbling mode) with the results obtained 
with the aid of the Marshalek moments of inertia (see 
Eqs.(1251),(IllD,(inni))- Evidently, if the "spurious" solution 
(the redundant mode) is not removed from Eq. (|19p . it 
contributes to the variables (|30p . These variables can not 
obey to the condition ([55)1 in this case. As a result, the or- 
thogonality condition Eq. (P7|) is broken and Eqs. ((5^ . (pni) 
for the transitions should produce completely different re- 
sults. A good agreement between the both calculations 
(see the right top panels in Figs lllll^ is the most valu- 
able proof of the validity of our approach. The observed 
negligible differences are due the approximate fulfillment 
of the conservation laws (PD]) . caused by the presence of 
the additional term ^ (see FigH]). 

According to our analysis, at hil ~ 0.25MeV there is a 
transition from the axially-deformed to nonaxial shapes 
with the negative 7-deformation in ^^^Yb (see Figl5]and 
the discussion in Refd. With a slight increase of the 




r = +1 r = -1 
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FIG. 12: (Color onhne) ^^^Dy. Similar to FigHU One 
observes a strong dominance of the B(E2)- and B(M1)- 
transitions from the wobbling states (r — —1) with spin I 
to the yrast states with spin I' — I ~ 1 starting from the 
rotational frequency hfl > O.SMeV. 



rotational frequency, starting from hft ~ 0.28MeV, the 
excited band of the negative signature, created by the 
first RPA solution, changes the decay properties. The 
negative values of 7-deformation produce the dominance 
of the inter-band quadrupole transitions from the one- 
phonon state to the yrast ones with a lower spin (A/ = 1, 
the case Eq. (|65|) . a)). 

Similar results are obtained in ^^^Dy for the lowest neg- 
ative signature excited band, created by the first RPA 
solution. At low angular momenta {Tifl < 0.3MeV) 
this band populates with approximately equal probabil- 
ities the yrast states with /' = / ± 1 (/ is the angu- 
lar momentum of the excited state). At Ti^l ~ O.SMeV 
a shape-phase transition occurs, that leads to the tri- 
axial shapes with the negative 7-deformation. In turn, 
the excited band, created by the first RPA solution, de- 
cays stronger on the yrast states with angular momenta 
I' — I — 1 (A/ = 1, the case Eg. ([55)1 . a)), starting from 
hfl > 0.32MeV. 

From the above analysis of the electric quadrupole 
transitions one notices that there is no need to know the 
definition of the wobbling phonon operator in the UR 
frame. Indeed, in this frame the direction of the angu- 
lar momentum is fixed and fluctuations of the angular 
momentum are absent. However, there is a vibrational 
mode, related to the shape fluctuations, that carries one 
unit of the angular momentum. According to Bohr and 
Mottelson consideration, in the PA frame the system 
shape is fixed, while the angular momentum fluctuates 
around the rotation axis that coincides with one of the 
principal axes of the inertia tensor. Evidently, the result 
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for the transition probabilities in the lab frame must be 
independent on the the choice of the reference frame. 

To prove the equivalence of the both results for the 
electric quadrupole transitions, let us use the Bohr- 
Mottelson definition of the wobbling phonon operator, 
Eq.dH]), 



)t _ 



2(4) 



(67) 

Here, the quantities W2.3 are determined by the Mar- 
shalek moments of inertia, Eqs. (P^ . In the PA frame 
one has to use the transformation Eq. (|A.12|l in order to 
calculate the transition probability and commutation re- 
lations Eq. (|A.2[) . This transformation and the definition 
(|67l) yield the expression 



B{E2;Iw 
1 



I±lyr) 



M. 



(E) 

2mi=± 



(68) 



that is similar to Eq. lpO)) . obtained in the UR frame, 
indeed. We will use this fact below to understand major 
features of the magnetic transitions from the wobbling 
band. 

In the CRPA approach the magnetic transitions are 
defined as 



BiMl-Jiy-^ I±lyr) 



(M) 
1Mi = 



(69) 



With the aid of the transformation from x-to z-axis quan- 
tization Eg. ([58)1 one obtains 



B{Ml;Iv 



/±l,r)«i|^[Afei,Ot]T[<ao,Ot 

(70) 

The linear boson term of the magnetic operator has a 
form (see also RefHT 



(71) 



where /xjv is the nucleon magnetons, gi'^^^\ g^'^^^^ are 
the spin and orbital effective gyromagnetic ratios, respec- 
tively; and the quasiparticle matrix elements aJT^^ are 
real. Taking into account the definition of the phonon op- 
erator Eq. (|15p . we express the magnetic transition with 
the aid of generalized coordinate and momentum ampli- 
tudes Eqs.m 



B{Ml-Jv^ I±lyr) 



(72) 



Are there any specific features of magnetic transi- 
tions from the wobbling band, related to the sign of 
7-deformation ? Our results evidently demonstrate the 
dominance of B{Ml; I nw I — lyr) (see the right bot- 
tom panels in Figs lllll2|) for the both nuclei. To get 
an insight of this result, we define the magnetic tran- 
sitions by dint of the Marshalek moments of inertia. 
With aid of the definition. Eg. ([57)1 . the transformation, 
Eg. (jA.12p . and commutation relations [Ii ± il2,M\m] = 
v/A(A + 1) - m(m T l)M\niTi for the PA frame [H], we 
obtain from Eg. (p5)) 



B{Ml-Iv I±lyr) 



(73) 



Although the expression for the magnetic transitions, 
Eg. ([75]) . is similar to the one of the Bohr-Mottelson 
model, we stress that the moments of inertia are defined 
self-consistently within the CRPA approach. Notice that 
the dipole magnetic moment (Af}*^^j^ [-1-]) increases quite 
drastically, if a nucleus is undergoing the backbending 
(see the discussion about Ml strength in Refl2l[). Keep- 
ing in mind that for the wobbling states W2.3 > 0, we 
have 

B{Ml;Inw ^ I-lyr) > B{Ml;Inw ^ I+lyr) (74) 

Thus, the tendency observed in the microscopic calcula- 
tions with the aid of the phonon amplitudes are under- 
stood in terms of the rules Eg. ([7^ . Independently on 
the sign of the 7-deformation of rotating nonaxial nuclei, 
these rules determine the dominance of A/ = ITi mag- 
netic transitions from the wobbling to the yrast states. 



IV. SUMMARY 

The observation of the wobbling excitations among ex- 
cited states near the yrast line is one of the most con- 
vincing manifestations of the nonaxial deformation in 
rotating nuclei. In this paper we present a transpar- 
ent, self-consistent derivation of the basic equations for 
the wobbling excitations in the UR (time-independent) 
frame, that determine: i) the energy spectrum and 
ii)electromagnetic properties of these states in even-even 
nuclei. We obtained the same expressions for the effective 
moments of inertia as those obtained by Marshalek 
in the time-dependent Hartree-Bogoliubov approach in 
the PA frame Q- We established one-to-one correspon- 
dence between the main characteristics of the wobbling 
excitations in the Bohr-Mottelson model and those that 
are derived within the CRPA approach. Note, however, 
that the CRPA breaks down at the transition point when 
Ap or A„ vanish [s^- We have avoided this problem by 
means of the phenomenological prescription for the rota- 
tional dependence of the pairing gap. A good agreement 
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between the dynamic moment of inertia calculated at the 
mean field approximation and the Thouless-Valatin mo- 
ment of inertia calculated in the RPA supports the con- 
sistency of our mean field calculations (see Fig[3]). In con- 
trast to the standard RPA calculations, where the resid- 
ual strength constants are fixed for all values of H. (see 
e.g. p^. Tiol. |20| ) ■ we determine the strength constants 
for each value of ^l by the requirement of the validity of 
the conservation laws. This enables us to overcome the 
instability of the RPA calculations at the transition re- 
gion, at least, for the excitations. In principle, projection 
methods may be used in the transition region in order to 
calculate transition matrix elements. Although the am- 
plitudes (j)'^^^ (see Eg. ([TCI) ) are higher for the RPA modes 
in the transition region than in other regions, the rela- 



tion < \'ipli''\ is stiU vahd to hold the QBA. And 

least but not last, the CRPA becomes quite effective at 
high spins, after the transition point, when the pairing 
correlations still persist. 

From our analysis it follows that an excited band can 
be considered as the wobbling one, if the magnetic transi- 
tions from this band into the yrast one fulfil the condition 
Eq. (|71)l . Note that these rules do not depend on the sign 
of 7-deformation. In contrast, the electric quadrupole 
transitions from this band to the yrast one must fulfil 
the staggering rules Ea. (p5)l depending on the sign of 
7-deformation. The transition from this band are char- 
acterised by a high collectivity. We predict that the low- 
est excited negative signature and positive parity band in 
^^^Yb (which is a natural prolongation of the odd angular 
momentum part of the 7-band) transforms to the wob- 
bling band at > O.SMeV. We found that strong E2 
transitions from this band populate the yrast states, with 
the branching ratio _B (£'2; /ly — > I —lyr)/ B{E2;I w 
I + lyr) > 1. Such a behaviour is brought about by the 
onset of nonaxial nuclear shapes with the negative sign of 
7-deformation after the backbending. According to our 
definition of 7-deformation, with the increase of the rota- 
tional frequency the system is driving to a noncoUective 
oblate rotation (around the x-axes). This trend is con- 
firmed by a good agreement between the results for the 
quadrupole transitions along the yrast line and available 
experimantal data (see Fig llOp . Similar transition occurs 
in i^^Dy after the backbending as weU, at > 0.3MeV. 
In this nucleus we found that the lowest negative signa- 
ture and positive parity band represents a natural pro- 
longation of the 7-band with odd spins, observed at small 
angular momenta up to Ml ^ O.SMeV. At rotational fre- 
quency Ml > O.SMeV this band created by the first RPA 
solution transforms to the wobbling excitations. A good 
agreement between calculated in the CRPA and experi- 
mental Routhians allows us to conclude that the exper- 
imental states, associated with (-I-, l)i band in Ref!27l, 
are wobbling excitations at the rotational frequency val- 
ues 0.45MeV< nn < 0.55MeV. These states fulfil all re- 
quirements that are specific for the wobbling excitations 
of nonaxially deformed rotating nuclei, with the negative 
7-deformation. Evidently, it is quite desirable to obtain 



new experimental data on electromagnetic decay prop- 
erties of these states in order to reach a final conclusion 
about our prediction and, respectively, on the validity of 
the CRPA analysis. 
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APPENDIX A: ASYMMETRIC ROTOR MODEL 

In this appendix we review the basic features of the 
wobbling excitations in the rotor model 1] in order to 
compare with the microscopic model discussed in section 
III. In addition to well known results, we provide a novel 
analysis of magnetic properties of the wobbling states. 

At high-spin limit / >> 1, the rotor Hamiltonian has 
the following form (see RefUl) 



23., (-^) PA ^" H^Qhi, 



H = Ho 

[I'i.,Ij]pA = -~iSijk{h)pA 



(A.l) 
(A.2) 



where (-/i)p^ and 5i are the angular momentum and 
principal moment of inertia components in the rotating, 
principal axis (PA) coordinate system. It is assumed that 
the yrast band is generated by rotation around the x-axis 
I m Ii = K . Small oscillations of the angular momentum 
create wobbling excitations described by the term 



(A.S) 



where 



W^2 = ^- ^, W3 = ^ - ^ (A.4) 

The wobbling (excited) state at spin / is created by the 
wobbling phonon 

OLbb = -yB ^ B+ = xd+^,,i + yO^obu (A.5) 
with 



(A.6) 
(A.7) 
(A.S) 



B+ = —{i2+ih)pA, B^{B^ 




(1^21^3) ' 
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and a normalization condition 



From 



Ea. (|A.3p it follows that the diagonalization of Hwobb re- 
quires W2 > 0, W3 > and 3i > ^2 > ^s- 

At {I ^ Ii >> 0) the eigenfunction of Hamiltonian 
(|A.1[) is a product of the Wigner V - function and the 
intrinsic eigenfunction defined by the wobbling quantum 
number = (Q^abb)"'"' l^)- ^he variable 

n„ = {IK\B+B\IK) K I-h=I -K (A.9) 

is defined with respect to the state \IK = {Ii)pa) = !"-«>) 
such that B+\IK) = ^/n^ + 1\IK - 1) = ^/nw + l\nw + 
1) and |0) = \IK = I). The transition probability for 
the operator of type X and multipolarity A 

B{XX; IvK^ I'v'K) ^ 

\{I'K'y\\M{XX)\\IK,i^)\^ (A.IO) 

is defined by the reduced matrix element 



{I'K'n'J\M{X\)\\IKn^) = VW+l x (A.ll) 
{n'J{IK\v^\I'K')M{X\v^ = /' - I)\n^) 



X 

E 

i/i — — A 



Notice that the eigenmodes of the Hamiltonian Hq do 
not change the projection K onto the first PA axis. How- 
ever, if the wobbling mode is excited with \n'^ — n^j] ^ 0, 
the projection K is changed and corresponding Clebsch- 
Gordan coefficients in Eq. (jA.ll[) can be expressed in 
terms of 5+, B operators (see Section 4.5c in Ref.l). 

Using this procedure for the Clebsch-Gordan coeffi- 
cients, a relation between multipole operators in the PA 
frame with the x- and z-quantization axis and the defini- 
tion Eq. (fTO| 



B{E2; In^^ I -2nyj) 



1 



(A.15) 

Here, we use the fact that (M{E2, vi = ±1)) =0 and 
{M{E2, vi = 2)) {M{E2, vi = -2)) and neglect terms 
of order or higher than ~ j. 

Similar procedure, described above, can be employed 
to derive Ml transitions from one-phonon wobbling band 
into the yrast band. Ignoring the terms of order I^^ in 
Eg. ljA.lip at high-spin hmit {I ~ K » 1), we have the 
following approximative values for the Clebsch-Gordan 
coefficients in terms of the matrix elements of the opera- 
tors B, B+ (or Qn,obb, Qtuobb) 



(I K 10 \ IK) 
{IK10\I±1K) 



1 

1 B+ 
±^(n±l|^ \n) 

T-^("T 1| ^+ \n) 



{IKl ±1\IK±1 
{I Kl ±1\I ±1K ±1) 1 (A.16) 

Since (7W(Afl t^i ±1)) = 0, we obtain 

{I'K'n'J\M{Ml,yi = 0)\\IKn^) = 

V27TT(nJ 6rj{M{Ml,iyi = 0)) + (A.17) 



+ Srj+i^{M{Ml,iyi = 0)) {xQt,M + yQ^oobb) 



M{XX^.,) = Y^ V^^^iO^l^O) M{XX,,l (A.12) jj^^^^ M{Ml-.,) = /..VSEti ( 



Yi=q\ii. 



we obtain for interband E2 transitions between the one- 
phonon wobbling band {n^ = 1) and the yrast band 
(n^ = 0) 

B{E2; In^-^I±lyr)^ T fcS^^^P (A.13) 

Here, the variables A^^^ and B^-^^ are defined by 



I i1=o]ii/i ) is a magnetic dipole op- 



erator. With the aid of Eqs. (TO]) . (rO|) . C02|) and the 

definition of operators given in Refi26, we obtain for Ml 
transitions from the one-phonon wobbling band {rim — ^) 
into the yrast band (riy, = 0) the following expression 



Eqs.lESD, (ESD, Qm ^ = e|Q™ and 



B{ Ml -In.^^ I±lyr) = 
1 (v^T\/W^)- 



= +!]> 



1 fw- 



1 fW'_ 



Q2 



/2I \W3j ' V2I \W2 

For intraband transitions we have 



(A. 14) 



47 yWWs 



(A.18) 



The magnetic moment {m['^^'Li[+]) can be calculated in 
any microscopic model. 
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